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Abstract 



f"^^ These notes present a critique of the standard three-flavor neutrino oscillation 

framwork. The design proposal of the MINOS at Fermilab based on a two mass 
eigenstate framework may require serious reconsideration if there is strong mixing 
between all three flavors of neutrinos. For the LSND and KARMEN neutrino os- 
dilation experiments, the amplitude of neutrino oscillation of the "one mass scale 
■ dominance" framework vanishes for certain values of mixing angles as a result of 

opposite signs of two equal and opposite contributions. Recent astronomical obser- 
vations leave open the possibility that one of the neutrino mass eigenstates may 
be non-relativistic in some instances. Neutrino oscillation phenomenology with a 
"p^^ superposition of two relativistic, and one non-relativistic, mass eigenstates is con- 

{2^ structed. It is concluded that if the transition from the non-relativistic to the rela- 

r-Gj tivistic regime happens for energies relevant to the Reactor and the LSND neutrino 

1/ . oscillation experiments then one must consider an ab intio analysis of the existing 

•^^^ data. 

.fe 

XI: 

1 Introduction 

I 

1 There are at least two aspects to neutrino oscillation phenomenology. The first 

is the kinematic aspect where neutrino oscillations are modeled in terms of 
mass squared differences associated with the underlying neutrino mass eigen- 
states and mixing angles. The exact form of the standard kinematic analysis 
depends on whether the weak neutrino eigenstates are of Dirac type, of Majo- 
rana type, or of the type recently proposed in Ref. [IH, and further on whether 
or not one allows for CP violation in the neutrino sector. The second aspect 
is dynamical in nature. This includes the well-known MSW modifications to 
neutrino oscillations [Q], the gravitationally induced neutrino oscillation phases 
H, and neutrino spin flips in strong gravitational and magnetic fields M. 
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d In this paper I shall confine my attention to the kinematical aspect of the 

^ phenomenon. The neutrino-oscillation parameters determined by kinematic 

^ analysis may then be used to determine which dynamical aspects may be 

Q important. It must be noted that even the kinematical aspects of neutrino 

^ phenomenology ^-0] are only now beginning to be understood properly 

1^ As soon as the preprint from Kamioka reporting the zenith-angle dependence 
of the atmospheric neutrino anomaly |Tl| arrived, it became clear that the 



O 

^ phenomenon of neutrino oscillations could no longer be studied in terms of the 

O two-fiavor analysis and that a great amount of confusion had arisen when the 

intuition and arguments from the two-flavor studies were carried over, without 

p I due reflection, to a realistic three-flavor analysis. Similar observations were 

O made by Learned, Pakvasa, and Weiler |Tp] (and few others) when they flrst 

^— > learned of the Kamiokande collaboration's initial results on the depletion in the 

. ^ flux of atmospheric muon neutrinos relative to the flux of atmospheric electron 

^ neutrinos. Due to conceptual and experimental complexities associated with 

D the neutrino-oscillation physics there are, at present, a large number of papers 

^ devoted to the three-flavor analysis. Still, some fundamental aspects of the 

j> three-flavor neutrino-oscillation analysis remain to be discussed. Some of the 

> remarks that I present are seemingly trivial, but in view of their possible 

relevance, I take the liberty of presenting them in this paper in the hope of 

Q establishing a better conceptual understanding. 

$-1 

Qh 

^ In the next section I briefly review the standard three-flavor neutrino oscilla- 

g tion framework. Section 3 is devoted to MINOS at Fermilab and other long- 

^ baseline experiments. Section 4 concentrates on the LSND neutrino oscillation 

experiment and KARMEN. Section 5 is conflned to a brief discussion of the Re- 
actor I7g experiments. Section 6 presents modiflcation to neutrino-oscillation 



00 

phenomenology with a superposition of relativistic and non-relativistic mass 
^sO eigenstates. □ Section 7 presents the conclusion. In addition to these seven 

sections, the paper contains three appendices. Appendix A provides further 

details on standard neutrino oscillation phenomenology. Appendix B contains 
j ) a simple theorem on the inverted mass hierarchy. Appendix C presents a ped- 

agogic discussion of neutrino oscillations and energy conservation. 

2 Three— Flavor Neutrino Oscillations: A Brief Review 



Because of its simplicity, the two-flavor neutrino oscillation framwork has 
dominated the design and understanding of neutrino oscillation experiments. 
However, as the observational and experimental indications for actual neutrino 



^ Section 6 is essentially based on a recent work of the author with Goldman 1 13 ] . 
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oscillations mount, newer and conceptually more sophisticated neutrino oscil- 
^ lation experiments must be considered. In these experiments, in my opinion, 

^ the insights gained from the two-flavor neutrino oscillation frameworks may 

Q prove to be inadequate. It is, therefore, important to have a fresh ab initio 

^ look at the three-flavor neutrino oscillations. I begin with a brief review of 

the three-flavor neutrino oscillation scenario and soon flnd myself embedded 
tkO in insights and results that would be either counterintuitive from a two-flavor 

^ neutrino oscillation point of view, or entirely absent unless the three-flavor 

O neutrino oscillations are invoked. Many of these results arise from a larger 

^ number of relative phases — relative phases that can lead, for example, to 

2 important cancellations. 

p I Assuming, (a) neutrinos to be Dirac particles, that (b) no CP violation oc- 

O curs in the neutrino sector, that (c) both the flavor eigenstates and the mass 

eigenstates are relativistic in the laboratory frame, and that (d) all three mass 
, , eigenstates are stable, the kinematically induced neutrino oscillation proba- 

^ bility, in a three-flavor neutrino oscillation scenario, reads: i 

> 



-I— > 



1 



VeAL,{r]k})= J dEf,iE)Pu'{E,L,{r]k}) , 



O 

^ where 

O 

$-1 

Oh 

\Q The kinematic phase that appears in the above equation is deflned as 

00 

^ a . L 



Pee' {E,L,{rik}) = See' - 4:Ue' lUeiUe' 2Ue2 sin^ (v^^i) 
-AUe'iUeiUe'aUessin^ (yl^) - AUe 2Ue2Ue' ^Ues sin^ [ifl^) ■ (2) 



V>l = ^^T7:z: ■ (3) 



\0 A 
I 

p/] Equations (|T]) and require several additional comments and observations. 

j ) These remarks are enumerated in Appendix A. However, the following deflni- 

^ tions need to be noted immediately: 

(i) The oscillations length, A^^^^, is deflned as 



The kinematic phase may also be written as: (p^^ = 1.27 Am^^ x (L/E). 
These kinematic phases may be modifled for dynamical reasons |@-^. 
Here, a = (5/2; (3 = 2.54 is the usual factor that arises from expressing 



'When no ambiguity is likely to arise I write Vee' {L, {%}) =V {ve ^ ^i' 
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d E in MeV, L in meters, and ^rri^jk • ^ refers to neutrino kinetic 

^ energy, a/p^ + m^, @ and L is the distance between the creation region and 

the detection region for the neutrino oscillation event. The five neutrino 



Oh 
=3 



O oscillation parameters {rfk] in equation (p are the two mass squared 

^ differences and the three mixing angles: ?7i = Am2i, 772 = Am|2, ''73 = ^, 

7]^ = (3, and t]^ = ip. The third mass squared difference is then given by 

^ (ii) fii^E) is the neutrino flux, of flavor £, at energy ii^. The fe{E) is normalized 

O to unity: 

D 

• 

O £ 

/ ciE/,(E) = l. (5) 

o 

-Emm and -Emax refer, respectively, to the minimum and maximum energy 
relevant to the neutrino beam and the neutrino detector. 
^ (iii) An element of the 3x3 unitary neutrino-mixing-matrix, U{9,P,ip), is 

^ labeled as Uij with i standing for any of the three neutrino flavors, i = 

^ e,fji,T] and j representing any of the three mass eigenstates, j = 1,2,3. 

The explicit expression for the mixing matrix that I use is (Maiani rep- 
resentation |T^, with CP phase 6 set equal to zero) 



O 

D 
> 

O 

$-1 

Oh 

u{e, (3, ^) 



ceC/3 secp sp 



,(6) 



in 

00 where cg = cos{6), sg = sin(6'), etc. 

1 

On 
I 

|3 3 The — > Ur MINOS at Fermilab, and Other Long-Baseline Ex- 

^ periments 

For the sake of concreteness let us choose a generic long-baseline z/^ — >■ Uj- 
experiment and consider the mass hierarchy such that 

'^3T2, 1) < ^ > -^2^ ^ L . (7) 



Or perhaps more precisely, E represents the expectation value of the kinematic 
Hamiltonian in the appropriate neutrino-flavor eigenstate. 
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Using Eq. (P and averaging sin^(- • •) associated with the ^3^(2,1) terms to 0.5 
^ I obtain 

Oh 

g - Vr) = Ac{P,ij) + Ao{9,P,i:) f dE U^{E) sin^ (^) , (8) 

where the ^-independent constant contribution, Ac{P, ifj), and the amphtude, 
^ Ao{9, P,ip), of the oscillatory term, in Vi^u^i Vt) are: 

O 

B 

O Ac{f5,i,) = 2 cos^(/3) cos2(V') sin2(V') , (9) 

^0(6', = 4 [cos(V') sin(6') + sin(/5) sin(?/^) cos( 
^-1 X [sin(-?/') sin(6') — sin(/?) cos(-?/') cos 

r-j X [sin(/3) cos('?/') sin(6') + sin('?/') cos 

X [sin(/9) sin(V') sin(6') - cos(V^) cos(6')] . (10) 

1^ Oftentimes one attempts to understand this experiment, and similar exper- 

P iments, within the framework of a two-flavor neutrino oscillation. Such an 

>■ analysis can be misleading. In the physically required three-flavor neutrino 

oscillations I see that there is a constant piece in V{Vfi Vr) (^nd an oscil- 
D latory term. Depending on the mixing angles the constant piece may make a 

O significant contribution to the V{v^ — > z/^) . 

Oh 

V{ue ^ yr) = BciP,^) + Boie,P,^) J dEUSE) sin^ {-f) , (11) 



The counterpart of Eq. (^) for the z/g Long-Baseline Experiments reads: 

'Co 



in 
\o 
00 

1 where 

P< i?c(/3,^) = 2 cos2(/3) sin2(/?) cos2(V;) , (12) 

Boie,P,ij) = 4 cos2(/3) cos(^) sin(^) 

X [sin('?/') cos(6') + sin(/?) cos('?/') sin(6')] 
hJ X [sin(V') sin(^) - sin(/3) cos(7/') cos(^)] . (13) 

Similarly, the counterpart of Eq. (|]) for the Ue — > i^fj, Long-Baseline Experi- 
ments reads: 

V{u, ^ u,) = Cc{P,ij) + Co{e,P,^p) J dEUSE) sin^ , (14) 

where 
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sin(?/') 



O 

o 

Cc(^,V') = 2 cos2(/5) sin2(/3) sin2(^) , (15) 

^ Co(^,/3,?/') = 4 cos2(/5) cos(^) sin(^) 

^ X [cos{%Ij) sin(6') + sin(/5) sin(-?/') cos(6')] 

2 X [cos(?/') cos{9) - sin(/?) sm(-?/') sm(^)] . (16) 

bX) It may be worthwhile to observe how, with respect to the mixing angle ^, a 

2 "large" ("small") 5c (/3,^) will complement a "small" ("large") Cc(/3,V)- In 

O particular, 6?/ measuring (in the same beam) the constant pieces in V^Ve ^t) 

■^-> and V{i'e — > t'/i) I can directly measure the mixing angles j3 and ip: 

> 

O V = tan-i , ^-^^^ , (17) 

5_i In Eq. ([18|) the angle is to be substituted from Eq. (|1 
> 
> 

O 

^ by measuring the amplitude of the oscillatory term in Vlu^^ — ^ z/^), or equiva- 

O lently any one of the Ao{9, P, ip), Bo{9, /?, ip), Co{0, (3, ip). Within the scenario 
represented by Eq. (0), the three probabilities (see comments below that re- 

'7^ quire an additional measurement) for i/^ — z/,-, z/e ^ z/,-, and z/g — z/^, measure 

^ four parameters — the three mixing angles [6, j3, ip) and one of the two in- 

1^ dependent mass square differences Arri^i. The long-baseline experiments, by 

\Q their very nature (within the scenario under consideration), determine only 
Am2i- The must then be determined either in a reactor experiment or 

1 an experiment like the LSND neutrino oscillation experiment. 

In principle, the three V{i'e — > z/£/) considered above can be measured at a 
very similar f^^{E). However, an additional measurement of one of the three 
h-p Vi^vi Vii) must be done at a sufficiently different fy^{E). This becomes 

^ obvious if I note that there are nine oscillation probabilities, 

P(z/^^Z/e), P(Z/^^Z/^), P(Z/^^Z/,), (19) 



Hence, the angle 6, the remaining of the three mixing angles, can be obtained 



that can be measured at a given energy E. Because of the three unitarity- 
imposed conditions, 

V{ui z^e) + v^) + V{yt ^ z/,) = 1, £ = e, /i, r , (20) 
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= V{u, - 





















1 



O 
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and the three conditions 

O 
Oh 

O 

C V[v^-^v;) = V[v,^v^) , (21) 

1^ arising from the tentative assumption that there is no CP violation in the 

"q neutrino sector, only three of the V{i'e //£/)[/(£')] are independent. There- 

^ fore, to determine the four parameters Am^i) one needs one more 

O measurement P(z/^ ^ ve)[f{E') ^ f{E)\. i 

> 

o 

4 Experiments Dedicated to — > z7e in the Appearance Mode: 
^ LSND and KARMEN 

3 



The published results of the LSND Neutrino Oscillation Experiment (NOE) 
and Karlsruhe-Rutherford Medium Energy Neutrino Experiment (KARMEN) 



both present results on v^. in the appearance (of i^e) mode |]I5|,|T6[. For 

b LSND NOE L ~ 30 m and for KARMEN L ~ 17 m. KARMEN, which is 

> less sensitive than LSND NOE, sees no neutrino-oscillation signal. The latest 

2 LSND NOE results, if interpreted within neutrino oscillation framework, yield 
Oh 

Xf^ For the /i"*" decay at rest, the normalized-to-unity Michel energy spectrum of 

^ the z7„ is given by 
00 

1 



O.SllE^;}^ ± 0.05 



X 10^1 



'^vi L 2 



^ My) = -^^-— , (22) 



I 



where y = 1/E, and 

7/2 = -J— ; = 52.8 MeV . (23) 



Substituting Eq. ( ^2] ) into Eq. (0) I obtain: 

Vi-^^^l^e) = - 4:U,iU^iU,2U^2 X 0{L,Aml,) 

- AUeiU^iUe3U^3 X 0{L, Ami,) 

- 4t/e2f/^2?7e3f/^3 X 0(L, Am^^) . (24) 



In Eqs. (^) to (21) all V{iy£ z/^) refer to same energy spectrum f{E). 
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=3 
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o 
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> 

a. 
o 

-I— > 
> 

> 

O 

D 
> 

O 

$-1 

Oh 

in 

00 
1 
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H-1 



In the above expression I have introduced a new function 0{L, Am^J. It is 
defined as: 



0{L, Aml) = 2yl 



C. 



+ 



c. 



sin (2 2/2 CJ^) + 



where C 



+ 2Cl^i{2y,C,;)- 



aL^m\ and 



Ci (2 2/2C,,) + 




Si(a;) 



sin(t) 



dt 



cos(t) 



(26) 



Given the mixing angles and detector efficiencies (that can depend not only 
on E but also on the location of an event within the detector), the function 
0{L, Am?) is a measure of the neutrino oscillation probability. I shall call 
0{L, Am?) the probability function. The probability function, corresponding 
to LSND's L ~ 30 m, is graphed in Fig. 1. 

As it should, 0{L, Am^) approaches 0.5 as Am^ oo. 

In an experiment using the neutrinos from the decay at rest of fi~^ (and in the 
absence of any directed magnetic fields), the event density within the detector 
is measured by the raw event density function, S{L, Am^), 



S{L, Am^ 



6(L) 

e(L) £\L, Am^ 



(27) 



where Am^ = {Am2i, Amgg} and e(L) is the energy-integrated average ef- 
ficiency at L inside the detector, and the definition of S'{L, Am^), the event 
density function, is obvious from Eq. (^). In actual experimental situations, 
such as LSND's NOE and KARMEN, the gain in event rate with decreasing 
L is inevitably shadowed by a increased background [|1^. By measuring the 
L-dependence of the event density function within a detector one may rule 
out one set of solutions over the other. 



To explore the physical content of Eq. (^) for Vl^f, 



let us consider 



degenerate [|^. Setting Am; 



an often-invoked scenario in which mass eigenstates \vxi and \v2i are almost 

- ^"^31 = ^"^3(2,1)5 and using Eq. (|), Eq. 

(28) 



2J) becomes: 

P(I^M ^ ^e) = A21 0{L, Aml^) + ^3(2,1) 0{L, Am2(2^,) 



with the two oscillation amplitudes defined as 
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^21=4 cos^(/?) cos(6') sm{9) cos(V') sin(6') + sm(/3) sin(V') cos(6') 



I 



I 

< 

H-1 



X 



cos('?/') cos(6') — sin(/3) sin('?/') sin(6') 



Oh 

O ^3(2,1) = 4 cos^(/3) sin^(/5) sin^(V') ■ (29) 

o 

In the scenario indicated above one usually assumes ^ m2 — mi. For 
'^"^21 ^ ^"^3(2 1)5 based on the Am^-dependence of sin^(t/CjJ, one 
O neglects the term associated with 1 and assumes that the dominant con- 

tribution to V{T'f^ —>■ Ve) comes from the Amg^-g^) term. This is the basic 
I scenario that the one-mass scale dominance framework considers. 

O 

I now make two observations in this context. 

D Observation I For (3 = mT/2 and/or ip = mc, n = 0,1,2, .. ., the oscillation 

'3 amplitude i .4.3(2,1) associated with the 0{L, Am^2 1)) term, the dominant 

^ term of the so called "one mass-scale dominance" identically vanishes. 

^— > Observation II Refer to Fig. 1, and consider a set of mass-squared differ- 

ences such that Am^, 1 lies roughly between 1 and 2 eV^ and Amg^g 1) is 

Q greater than 2 eV^. Then 

$-1 

Oh 

in 

^ For the above configuration of mass-squared differences, Eq. ( |30D violates 

^ one of the assumptions of the "one mass-scale dominance." That is, the > 

Q\ in the above expression is in contradiction to the expected <^ of the the one 



0{L,Aml,) > Am2(2,i)) . (30) 



mass-i 



scale dominance analysis. @ 



^For "and" in "and/or" above the three-flavor analysis reduces to a two-flavor 
framework. 

^ Note should be taken that specific values of Am^ considered here refer to L = 30 
meters for LSND NOE. In addition, as the detector has a width of about 8 meters 
centered at L = 30 meters, depending on the oscillation length under consider- 
ation, appropriate L-integration of Eq. (^) may need to be considered. Similar 
considerations apply to KARMEN. 
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Observations I and II suggest that the 0{L, Arn^i) term may make a domi- 



nant contribution to V{i'^ 



for certain values of P, and Am; 



2 

21- 



< 



0.8 - 



0.6 



0.4 ^ 



0.2 ^ 




Am^ (eV^) 

Figure 1. The function 0{L, Am^), 
LSND's Lc^30 m. 



corresponding to 



For the indicated-example configuration of mass-squared differences, the canon- 
ical wisdom of the one-mass scale dominance fails, in part, because V{l7fj, —>■ 
Ve) evaluated at 



average 



3/2 



hSy) 



y 



(31) 



does not equal Michel-spectrum averaged Viv^^ — > z7e). This is particularly 
true for low Am^, that is, before the Am^-region where 0{L, Am^) reaches 
its (Am^ —>■ oo)-value of one half [which in turn coincides with the average 
value of sin^(- ■ ■))]. 

The above discussion should not be interpreted to mean that the existing 
neutrino-oscillation data necessarily imply the specific values of parameters 
for which the one-mass scale dominance requires revision. It remains possible 
that actual neutrino-oscillation parameters satisfy the canonical one-mass 
scale analysis. In fact, the LSND NOE relevant Am^ of Ref. [1^ lies within the 
applicability of one-mass scale dominance analysis as will be shown elsewhere. 
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However, it is important to define the unexpected boundaries where the one- 
mass scale dominance analysis can, and does, break down. 



e 

S 0.8 

B 

< 0.6 
(J 

b 



0.4 - 



0.2 - 



L=30.0m, E^j„=20.0 MeV 




10 



Am^ (eV^) 



10 



Figure 2. The functions C(L, Am^ , E^in) for LSND 
NOE. 

For various reasons, such as detector efficiency, neutrino background, etc. 
|T3|,|TB|, one may wish to introduce an energy cutoff in E^^. For LSND's decay- 
at-rest data this energy cutoff corresponds to Emin — 20 MeV. Under such 
circumstances, assuming that there is no cut on E^ax except that dictated by 
the Michel spectrum itself, the 0{L, Am|J in the above discussion should be 
replaced by 



0\L, Am^ , E 



-1 r 



4^/2 



y2 



2yi 2yt 



1 

2/2 



+ sin (2 y2 C^J - 



yi 



sin (2 yi C^J 



H cos (2 y2 Cj,) 

1/2 



cos (2 yi CjJ - 7^ cos (2 ?/i C^J 



2 2/1 



+ 3 C,, Si (2 y2 C,,) - 3 C,, Si (2 C,, 



+ 2 ^2 C^,'. Ci (2 y, - 2 Cl Ci (2 C^,, 



(32) 
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where yi = 1/Emin- The new function 0'{L, Am^^, Emin) is plotted in Fig. 2 
for LSND. The differences and similarities between the functions 0{L, Am^^, 
and 0'{L, Am^J are apparent in Fig. 2. 



5 Experiments Dedicated to z/g Ve in the Disappearance Mode: 
Reactor V^, Experiments 



These experiments are performed with reactor I7g, which have a typical energy 
of about 5 MeV. For this section I shall assume E = 5 MeV. In extracting the 
five neutrino-oscillations parameters from the existing data we must use the 
exact known spectra and detector efficiencies. 

At the present time there exist at least four results on Vl^e ^e) from reactor 



experiments |l9|-p^. These experiments work in the Ve disappearance mode. 
The core and the reactor sizes place a lower limit on the distance between 
the core and the detector. For the latest experiment from Savannah River 
we have L of about 18 m, for one position of the detector, and roughly 24 m 
for another. The overall best limit on VlVf. —>■ z7g) comes from the Bugey 



experiment pT[. Results from all reactor experiments are consistent with the 



null-oscillation hypothesis. The oscillation hypothesis, within the framework 
of a two-flavor neutrino oscillation, for the Savannah River experiment |2^ 
yields Am? = 3.84 eV^. The existing cosmological arguments suggest that 
Hif^i < 30 eV. Thus, the range eV^ < Am^^ < 10^ eV^ holds particular 
interest. 

There are two independent Am^^. Let these be Am\^ and Amgg. Then, the 
third Am^j is simply given by: Am\^ = Am^i + Am|2- Associated with these 
Am^j are three oscillation lengths (of which only two are independent): Ag^^, 
, and 



V osc 



\osc \0 

A-= ■ (33) 

\OSC I \osc ^ 
^21 + ^3 2 

The initial three independent length scales associated with the masses of 
the mass eigenstates reduce to two independent length scales. The physical 
origin of this fact lies in the assumption that all three mass eigenstates are 
relativistic. 

To gain some understanding of the neutrino-oscillation parameter space 
(Am2i,Am32, 6', /?,■?/') suppose that Am^i ~ Am32 ^ 10^ eV^. Then, A2T ^ 
Agf ~ 3 cm and a detector placed at L 3> 3 cm sees only an overall z/g deficit: 

V{Ve ^ I7e) = 1 - 2 cos^{(3) (sm\f3) + cos2(/?) cos2(^) sm\e)) . (34) 
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This expression follows by setting all sin^(- ■ ■) in Eq. (|T]) equal to 0.5, as is 
appropriate when the corresponding A°'''^ ^ the relevant L, setting £ = = e, 
and exploiting certain trigonometric identities. 

For this scenario, the reactor experiment is unable to determine any of the 
Am^j. It is only sensitive to the two of the three mixing angles. The z/g deficit, 
if Am2i — Am|2 — 10^ eV^, in reactor experiments is independent of the 
mixing angle ip. The Savannah River experiment yields (in the "integrated 
rate test") P(z7g v^) = [98.7 ± 0.6 (stat.) ± 3.7 (syst.)] x 10"^ for "position 
1" of the detector and Viu^ V^) = [105.5 ± 1.0 (stat.) ± 3.7 (syst.)] x 10"^ 
for "position 2" of the detector. These rates are in good agreement with the 
no-oscillation hypothesis. 

Since, in the above defined scenario, the condition L ^ 3 cm for the Sun is 
satisfied by thirteen orders of magnitude, and if CP is not violated in the 
neutrino sector, and a solar neutrino deficit arises solely from a kinematically 
induced neutrino oscillation, then Eq. (|34D is valid for Sun as well. That is, in 
the above defined scenario, P(Fe ^ i^e)|Reactor = 'Pi'^e '^e)|soiar- But since 
|23| P(z/e — ^ t'e)|soiar ~ 0.5, the sceuario Am^i ^ Am|2 — 10^ eV^ is ruled 



out. So, at least one of the Am^^ <^ 10 eV . 

This is a generally known conclusion but, to the best of our knowledge, here 
it has been reached via an entirely new chain of arguments. 



6 Neutrino Oscillations with a Superposition of Relativistic and 
Non-Relativistic Mass Eigenstates 



The existing indications of neutrino oscillations [^,0,^ arise from the data 
that contains neutrino energies in sub-MeV to GeV range. The neutrino os- 
cillation phenomenology within the standard three-flavor framework contains 
the fundamental assumption that neutrino mass eigenstates that superimpose 
to yield the neutrino flavor eigenstates are relativistic. This assumption is sup- 



ported by the cosmological argument that the sum of the neutrino masses 
have an upper bound of about 30 eV 



J2 H^^e) < 30 eV 



(35) 



However, very recent astronomical observations ||25|-p7| raise potentially seri- 
ous questions on the validity of the standard cosmological model. First 
the UC Berkeley's Extreme Ultraviolet Explorer satellite's observations of the 
Coma cluster of galaxies indicates that this cluster of galaxies may contain a 
submegakelvin cloud of ~ lO^^M© baryonic gas. Second [^, the discovery by 
the German x-ray satellite Rosat in which a sample of 24 Seyfert galaxies con- 
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O 

o 

tained 12 that were accompanied by a pairs of x-ray sources, almost certain to 
^ be high redshift quasars, ahgned on either side of the galaxy. These observa- 

^ tions and the associated interpretations, if correct, may place severe questions 

Q to cosmological models that depend on the ratio of photon to baryonic density 

^ in the universe and its size. 

Tentatively, therefore, I relax the cosmological constraint (|35|) completely and 
1-) explore the resulting consequences for three-flavor neutrino osccillation frame- 

"q work. The latest kinematic limits on neutrino masses are much less severe 

O 

m{vr) < 23 MeV , (36) 

m(z/^) < 0.17 MeV , (37) 

-S m{ue) < 10 - 15 eV . (38) 

^ The existing and the proposed neutrino oscillation experiments involve neu- 

^ trino energies from a fraction of a MeV, if not less, to several hundred GeV for 

the upper-energy end of the neutrino beams. It remains possible that for some 
5-1 of the experiments (or a certain sector of an experiment) the mass eigenstates 

^ lz/3) and lz/2) are non relativistic. Towards the end of understanding such a pos- 

O sible situation I now consider the interplay of non-relativistic and relativistic 

mass eigenstates in three-flavor neutrino oscillation framework. 

O 

$-1 

Qh Consider a physical situation where we have two relativistic, \vi) and lz/2), and 

one non-relativistic, jz/a), neutrino mass eigenstates. 

At t = 0, X = 0, assume that a source creates a vp 

in 

00 k^) = Un\i^i) + Ut2\y2) + t^aks) • (39) 

The spatial envelope, which is assumed to be "sufficiently" narrow, associated 
<^ with {|z/i), 11^2)} evolves towards the detector as a; ~ t, while the spatial 

|V| envelope of evolves towards the detector as 



< X ~ (p/ma) t = t . (40) 

Therefore, the {|z/i), |z/2)} arrives at the detector at time tj ^ L, while the 
1^3) arrives at the detector at a time tjj = m^L/ ^2m3{E — m^). For the 
above considered energies and for a source-detector distance of a few tens (or 
greater) of meters the detected has no overlap with at the registered 
event at t/, and similarly the detected |z/£/) has no overlap with {jz^i), 1^2)} 
at the registered event at t//. With these observations at hand one can eas- 
ily evaluate the modification to the neutrino oscillation probability (H). The 
modified expression reads: 
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o 
o 

g Pee {E, L, {^,}) = {U^z U^^f + [ [Uei + U.^f 

^ - AU,,^UnU,2U,2^i^H'^lM . (41) 

5_( The first term on the rhs of the above equation is the contribution from the 



I 

< 

H-1 



First: 



^ non-relativistic mass eigenstate to the Pu' {E, L, {^k}) at the ue' event at 

•gj^ time t// while the second term is the contribution from the relativistic mass 

(jj eigenstates to the i>£i event at the earher time tj. Contained in the above 

expression is the the fundamental "collapse of the wave packet" postulate of 

D the orthodox interpretation of the quantum mechanics. That is, given a single 

Q z/£ emitted at the source, if the event occurs at tj no event occurs at tu, and 

$—1 

p ^ vice versa. 

> 

p I 

Q There are several observations that one may make about the result (pl). These 

^ observations follow. 

The neutrino oscillation probability now contains only one length scale, 
5_( Ci = ^^^21 5 ^2 = 6', ^3 = /?, and ^4 = However, this loss of length 

^ scale is related to a manifestly different expression for neutrino oscillation 

O probabilities. 

D 
> 

P Second: 
Oh 

If 777-3 is in the range of a fraction of an MeV to a few tens of MeV, one 
^ cannot base the analysis of existing neutrino oscillation data in terms of Eq. 

1^ (0), or Eq. (^TJ), alone. For the zenith-angle dependence of the atmospheric 

\Q neutrino anomaly data |ll|] neutrino energies in the GeV range are involved. 

This meets the requirements under which Eq. (Q) is operative. On the other 
hand, the physical situation for part of the LSND events (energy range 
^ between 20 MeV and Michel spectrum cutoff of 52.8 Mev JT^), and the 

I reactor experiments (average z/g energy about 5 MeV 0-|22|) (^I)) , the 

physical conditions for Eq. (¥T) may be satisfied. 



Third: 

For the solar-neutrino deficit one may speculate that 777.3 niay be such that 
an energy-dependent transition takes place from the non-relativistic regime 
to relativistic regime, thus accounting for the apparent energy dependence 
of the solar-neutrino deficit. So, consider a situation where the length scales 
are such that all sin^ iy^i) average to 1/2. Then 

Pee[Eq. (I)] = 1 - 2U!,U!, - 2Ul,Ul - 2U!,U!, , 

Pee [Eq. (Ell)] = f/,\ + + t^e3 ■ 
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o 

d Exploiting unitarity of the mixing matrix 1 immediately see that above 

^ speculation, within the defined context, has no consequence because 

Oh 

O 

^ In astrophysical environments, such as for neutrinos observed in the super- 



Unitarity P^e [Eq. (§] = P^e [Eq. 
Fourth: 



> nova 1987a |^O|,0, it may happen that is such that the evolution of 

^ envelope does not escape the astrophysical environment. That is, p/m 

associated with the envelope of It's) is less than the escape velocity 

P-i r 

^ Non Relativistic: E < m + — , (42) 

M-i 2 r 

O 

where Tg = 2GM is the gravitational radius of the astrophysical object. 
]^ Under these circumstances Eq. (^Tj), for a detector at Earth, reduces to: 

^ Pw {E, L, {a}) = e (^E - m + ^) 

2 + + - 4 f/,, 1 f/£ 1 f/^' 2 f/£2 sin^ ((^0 J . (43) 

O In the above expression 0(- ■ ■) is the usual step function, vanishing for its 

-t— > argument less than zero and equal to unity for its argument greater or equal 

2 to unity. 

Oh 

Now 1 rewrite Eq. ( ^T]) in a form that makes the deviations of result (|41|) 
^ explicit from the corresponding two-fiavor scenario with two relativistic mass 



ly-^ eigenstates. This new form of Eq. (|4T|) reads: 

^ Pw {E, L, {a}) = - sin^ ((^0 J , (44) 

1 

2 2 

with A^^/ = {Ui'3 Ues) + (t/^'i t/^i + t/F2 ^^£2) , or more exphcitly 

I 
I 

< 

H-1 





^^ + c^-l 


Op2 2 2 


9p2„2 2 


\ 


A = (5 + 


2p2 2 2 


2c^s^(c^s^ — 1) 








Op2 2 2 

y ^c^s^c^ 


Op4 2 2 


2 2 / 2 2 _ 





(45) 



and 

=4f/^,lt/£lf/£'2^£2 ■ (46) 

In Eq. (P5|) 5 is a 3 x 3 identity matrix. 



Only when both (3 and ■?/' vanish does Eq. (^4]) reduce to the expression for the 
two flavor scenario with two relativistic mass eigenstates — for then U{9, /3, ip) 



16 



o 
o 

o 

Oh 

=3 

O 

o 
o 

• 

o 

> 

a. 
o 

-I— > 
> 

> 

O 

D 
> 

O 

$-1 

Oh 

in 

00 
1 

On 
I 

I 

< 

H-1 



becomes block diagonal with no mixing with ji/s), as 



^ 1 0^ 



p=^=0 : A 



1 
1 



/ 



A 



"^2 
^26 



^29 U 



^26* ^26 



V 











(47) 



It is, therefore, concluded that there is enough richness in the three flavor neu- 
trino oscillation phenomenology that unless an argument can be made that 
this structure is incapable of accommodating the existing neutrino oscillation 
data there is no necessity to put forward the sterile neutrino hypothesis. I 
hasten to add that the purpose of this paper was not to show that the un- 
covered structure and length scales can necessarily accommodate the existing 
and the forthcoming data, but only to reveal the often ignored structure and 
the hidden complexity of the three-flavor neutrino oscillation phenomenology. 



7 Conclusions 



In these notes I presented a critique of the standard three-flavor neutrino 
oscillation framwork. I have argued that the standard three-flavor neutrino 
oscillation framework has rich structure that must not be ignored in making 
various approximations and designing experiments. Specifically, the design 
proposal of the MINOS at Fermilab based on a two-mass eigenstate framework 
may require serious revision if there is strong mixing between all three flavors 
of neutrinos. The general lesson to learn here is that: 

Even if one is searching for neutrino oscillations between two flavors it may 
not be advisable to work in terms of the well known formula (cf, [0, p. 17]): 



Pii^a i^h) = sin2(20) sin^ 1.27 



E 



where Am is measured in eV , L in meters (kilometers), and E in MeV (GeV). 

Detailed analytical and graphical analysis of LSND NOE revealed the interest- 
ing roles that various relative phases play in neutrino oscillation probabilities. 
For example, in the context of LSND NOE and KARMEN, the amplitude 
-^3(2,1) 5 defined in Eq. (pQ]), vanishes for certain values of mixing angles as a 
result of opposite signs of two equal and opposite contributions, i I further 



Such cancellations are a generic feature of three, or more, state systems and range 
from neutrino oscillation phenomenology to laser oscillation without population 
inversion [bSI. 
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argued, for instance, why the canonical wisdom of "one mass scale dominance" 
^ may fail in certain unexpected sectors of the neutrino-oscillation parameter 

^ space. By comparing the solar neutrino deficit and the results on 17^ disappear- 

Q ance in the reactor experiments I reached the conclusion that the standard 

^ three-flavor neutrino oscillation framework requires one of the mass-squared 

differences to be -C 10^ ev^. This, when coupled with the very recent astro- 
bX) nomical observations, left open the possibility that one of the mass eigenstates 

may be non-relativistic and hence introduce a fundamental modification to 
O the standard neutrino oscillation phenomenology without invoking a sterile 

^ neutrino. If the transition from the non-relativistic to the relativistic regime 

2 happens for energies relevant to the Reactor and LSND NOE then one must 

_Ph consider an ab intio analysis of the existing data. 

Oh 

O Acknowledgements 

I extend my thanks to Drs. Christoph Burgrad, Terry Goldman, Peter Herczeg, 
^ Mikkel Johnson, Bill Louis, John McClelland, Ion Stancu, Hywel White, and 

Nu Xu for our continuing conversations on neutrino oscillations and other 
g matters of physics. 

^ This work was done, in part, under the auspices of the U.S. Department of 

O Energy. 

O 
$-1 

Oh 

in 

00 

I 
I 

< 



18 



A Comments and Observations on Eqs. (1) and (2) and A Brief 
Outline of the Standard Neutrino— Oscillation Phenomenology 



Equations (|l|) and (H) require several comments and observations. These re- 
marks follow. 

While analyzing experimental situations, the integral appearing in Eq. (^) 
is evaluated by summing over approximately-constant bins of fe{E). In this 
context 1 note that the integral 

h= J rfEsinM2vr^ = / rfE sin^ (^) , (A.l) 

with Cjt = aLAm^j may be evaluated by introducing y = 1/E, and using the 
result 



dE sin^(%)=/(-Qsin^(,C, 



cos(2,,C,.) _ 

2 2/2 ^ " ' iy2 



+ ^^^^%^ + Si(2,,C,)C,--^ . (A.2) 
^ y 1 ^ y 1 



Referring to Eq. @, it is important to note how it arises. Apart from the 
already indicated assumptions, the fundamental assumption of the neutrino- 
oscillation phenomenology, contained in Eq. (^, is that weak flavor eigenstates 
for i/g, v^^ z/^, and I7g,l7^,!7^, are not mass eigenstates. Why this is so, no one 
knows at the present time and the question remains a deep mystery. To ap- 
preciate the nature of this mystery, note that from the 1939 classic paper of 



Wigner we have learned that every physical state is described by two 
Casimir invariants associated with the Poincare group of space-time sym- 
metries. These Casimir invariants are directly connected to mass and s'pin. 
Why, therefore, do weak interactions and space-time symmetries intermingle 
in such a manner as to prefer linear superpositions of mass eigenstates? I shall 
return to this subject elsewhere building on some recent work on the subject 
|]J3|^. 



Here is a brief outline of the standard neutrino-oscillation phenomenology. 
Assume that in the "creation region," Tic, located at Vc, a weak eigenstate with 
energy E (with appropriate uncertainty dictated not only by the Heisenberg 
Principle but also determined by the masses associated with mass eigenstates) 
denoted by TZc), is produced at Vc, with the clock set to t = 0. Each of the 
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three neutrino mass eigenstates shall be represented by |z/j); i = 1,2,3. So I 
have the linear superposition: 

We, T^c) = J2 W^) ' (^-3) 

1=1,2,3 

where i = e, fi,T represent the weak flavor eigenstates (corresponding to elec- 
tron, muon, and tau neutrinos, respectively). The {ui), |z/2), and lz/3) corre- 
spond to the three mass eigenstates of masses mi, m2, and m^, respectively. 
Under the already-indicated assumptions the unitary mixing matrix Ug^ may 
be parameterized by three angles and is given by Eq. (|^). 

At a time "t"= t > 0, I wish to study the weak flavor eigenstate in the 
"detector region," TZd, located at ra- Under the above indicated assumptions 
the neutrino evolution is given by the expression 



|7^rf) = exp 



-'-jmi+'-jp-dx\\u,,n,) . (A.4) 



Here H is the time translation operator, the Hamiltonian, associated with the 

— * 

system; P is the operator for spatial translations, the momentum operator, 
and H(t,x), P(t,x) = 0. It shall be noted that |z/^. Tic) has evolved to a 
state \7ld) that, in general, is not an eigenstate associated with z/g, u^, or z/,-. 
Instead, it is a hnear superposition of states associated with z/g, u^, or z/,-. 
The "neutrino oscillation probability" from a state \i>i, TZc) to another state 
We', Tld) is now obtained by calculating the projection (z/^/, lZd[R-d), i-e., the 
amplitude for |z/^, TZ^ \i>£', TZd), and then multiplying it by its complex 
conjugate. An algebraic exercise that exploits (i) the unitarity of the neutrino 
mixing matrix U{6, f3, (ii) orthonormality of the mass eigenstates, (iii) 
certain trigonometric identities, makes (iv) the standard observation (for neu- 
trino oscillations) that {Ej — E^) = Am^^c^/ (2^^), where Am^^ = m| — and 
~ = E, and (v) sets \rd - = L, yields Eq. (g) for P«/ {E, L, {r/^}). 



B A Simple Theorem on the Inverted Mass Hierarchy 
Under 

A^2l"A^32 , (B.l) 

and simultaneous change in U{6, f3, tp) of the form 

f/^'if/n "±f/^'3t/,3 , (B.2) 

Ue>2Ue2^ — >±Ui'2Ui2 , (B.3) 
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o 

o □ 

the V{iy£ — > ugi), given by Eq. (P, remains unchanged, u 

O 

Ph The invariance of various neutrino oscillation probabilities Vl^e ue), under 

Q the above defined change in neutrino-oscillation parameters, does not imply 

^ that other physical phenomena too remain unaffected under the change in the 

^ neutrino-oscillation parameters defined above. 

O For considerations on the evolution of the universe, and formation of the large- 

"q scale structures in the universe, it may be noted that given a mass eigenstate 

^ configuration ^ 1712 > mi (by definition, the "standard mass hierarchy"), 

O the corresponding "inverted mass hierarchy" (by definition, > m2 S> mi) 

p I defined by the above symmetry, necessarily has a larger than the 

"p_i associated with the (corresponding) standard mass hierarchy. 
O 

Within the above framework (i.e., three fiavor neutrino oscillation phenomenol- 

J3 ogy with Dirac neutrinos without CP violation), the neutrino-oscillation ex- 

^ periments cannot distinguish between the standard mass hierarchy and the 

(]_) inverted mass hierarchy. Additional physical phenomenon must be considered 

^ to distinguish the standard and inverted mass hierarchies. 

> 



C Energy Conservation and Neutrino Oscillations 



O 

D 
> 

O 

$-1 

Oh 

This section is included here purely for pedagogic reasons, and may be skipped 
^ by the "experts." 

in 

\0 Let us begin with posing a pedagogic question. To keep the algebraic details 

^ simple, I will confine to a situation where the mixing angles (3 and ip vanish. 

1 Then, consider an electron neutrino in the creation region: 

I 

and, again for the sake of argument, I will confine to a situation where the 
neutrino mass eigenstates have same momenta, and hence different energies. 
The expectation of value of energy in the state given by (|C.1|) is 



|z/e, 7^c) = ce\vi) + se\v2) , (C.l) 



{E,:)=clEi + slE2 . (C.2) 



^ Either the upper + sign, or the lower — sign, in both of the above equations is 
meant to be taken in any given transformation. 

Similar comments apply in the context of gravitational effects on physical phe- 
nomenon involving neutrinos in astrophysical environments. An exception is the 
gravitationally induced modification to the neutrino-oscillation probabilities, where 
a similar symmetry as mentioned above exists ||3|. 
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The posed question is: Where did the excess energy, 



O 

o 

d The state |z/e, TZc) evolves with time, and at any given time it has a definite 

^ probabihty of being found (on measurement) as an electron neutrino, and 

^ one minus that probability of being found as a muon neutrino. The energy 

Q expectation value of the energy remains constant in time. But now suppose I 

^ make a measurement and detect a muon neutrino. The linear superposition of 

Ch the electron and muon neutrino state now collapses to muon neutrino 

^ W,., = -se\ui) + ce\u2) . (C.3) 

> 

O 

The expectation value of the energy is now 
2 {E,^)=slE, + clE, . (C.4) 

O 

^ AEexcess = {EuJ - {E,) = C2e iE2 - E,) , (C.5) 

^ come from? 

^ To answer this question it should be realized that in order for the detector to 

^ detect a muon neutrino it must be able to make a measurement over a period 

t 

^ At ^ —V = . . , ^ = ^^1 + ^2^2 . (C.6) 

Oh 

As such, the act of measurement must impart a minimum energy (via the 
^ scattering process in the detector) AE Heisenberg =~ fi/{2At) to the detected 

. ^ muon neutrino. It is readily found that 

in 

'-^-L^ Heisenberg — ■ {^■' ) 

0\ Referring to Eq. ( |C3D , now note that E2 - Ei ^ Aml^c^/{2 E). Thus 

I 

So what at first might have appeared as a violation of energy conservation is 
precisely the energy, or more accurately the C2e fraction of AE Heisenberg, that 
the act of measurement imparts to the detected muon neutrino. 



AEexcess — C29 AEneisenberg ■ (C'-S) 



For a more formal discussion of this point see Ref. 
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In this section I shall not set h and c equal to unity. 



22 



References 



[1] D. V. Ahluwalia, Int. J. Mod. Phys. A, 11, 1855 (1996). 

[2] S. P. Mikheyev and A. Yu. Smirnov, Sov. J. Nucl. Phys. 42, 913 (1985); 
L. Wolfenstein, Phys. Rev. D 17, 2369 (1978). 

[3] D. V. Ahluwalia and C. Burgard, Gen. Rel. and Grav. 28, 1161 (1996). A 
continuing debate exists on the subject of "gravitationally induced neutrino 
oscillation phases." In our opinion, our reply and the preprint of Y. Grossman 
and H. Lipkin is sufficient to resolve the confusion. The relevant references are: 
T. Bhattacharya, S. Habib and E. Mottola, |gr-qc/960507^ ; D . V. Ahluwaha 
and C. Burgard, gr-qc/9606031 ; Y. Grossman and H. Lipkin, ^r-qc/9607201 



C. Y. Cardall and G. M. Fuller, |hep-ph/961049^ ; N. Fornengo, C. Giunti, C. 
W. Kim, and J. Song, [hep-ph/9611231 . The basic result of our work can be 
summarized as follows: The phenomenon of neutrino oscillations provides a 
"flavor oscillation clock," and this clock red shifts as required by Einstein's 
general relativity. 

[4] D. Pfriz, M. Roy, and J. Wudka, Phys. Rev. D 54, 2761 (1996). 

[5] C. M. Kim and A. Pevsner, Neutrinos in Physics and Astrophysics (Harwood 
Academic Publishers, 1993). 

[6] B. Kayser (with F. Gibrat-Debu and F. Perrier), The Physics of Massive 
Neutrinos, (World Scientific, Singapore, 1989). 

[7] R. N. Mohapatra and P. B. Pal, Massive Neutrinos in Physics and Astrophysics 
(World Scientific, Singapore, 1991). 

[8] H. J. Lipkin, Phys. Lett. B 348, 604 (1995). 



[9] T. Goldman, \hep-ph/9604357 . 

[10] E. Alfinito, M. Blasone, A. lorio, and G. Vitiello, Phys. Lett. B 362, 91 (1995); 
Also see, LANL preprint archive: hep-ph /9601354 . 

[11] Y. Fukuda et al., Phys. Lett. B 335,237 (1994); 

For a comment on the statistical significance of the quoted zenith-angle 
dependence of the atmospheric neutrino anomaly by Fukuda et al., see D. 
Saltzberg, Phys. Lett. B 355, 499 (1995); 
B. C. Barish, Nucl. Phys. B (Proc. Suppl.) 38, 343 (1995). 

[12] J. G. Learned, S. Pakvasa, and T. J. Weiler, Phys. Lett. B 207, 79 (1988). 

[13] D. V. Ahluwalia and T. Goldman, forthcoming. 

[14] L. Maiani, Phys. Lett. B 62, 183 (1976). 

However, the manner in which rows and columns are labeled and notation for 
angles is essentially that of Ref. ||5|, Eq. 6.21]. 



23 



[15] C. Athanassopoulos et al., Phys. Rev. Lett. 75, 2650 (1995); 

C. Athanassopoulos et al., Phys. Rev. C (1996, submitted); LANL archive 
preprint: flucl-ex/9605001 . 



[16] B. Bodmann et al, Phys. Lett. B 267, 321 (1991); B. Bodmann et al, Phys. 
Lett. B 280, 198 (1992); B. Zeitnitz et al, Prog. Part. Nucl Phys. 32, 351 
(1994). 

[17] W. C. Louis, (private communication, April 1996). 

[18] G. L. Fogh, E. Lisi, and G. Scioscia, Phys. Rev. D 52, 5334 (1995); 
G. L. Fogh and E. Lisi, Phys. Rev. D 52, 2775 (1995). 

[19] G. Zacek et al., Phys. Rev. D 34, 2621 (1986). 

[20] G. S. Vidyakin et al., Pis'raa Zh. Eksp. Teor. Fiz. 59, 364 (1994) [Enghsh 
translation: JETP Lett. 59, 390 (1994)]. 

[21] B. Achkar et al., Nucl Phys. B 434, 503 (1995). 

[22] Z. D. Greenwood et al., Phys. Rev. D 53, 6054 (1996). 

[23] Results from SAGE: J. N. Abdurashitov et al., Phys. Lett. B 328, 234 (1994); 
Results from GALLEX: P. Anselmann, Phys. Lett. B 327, 377 (1994); 
Results on ^B Solar Neutrinos from Kamiokande II: K. S. Hirata et al., Phys. 
Rev. D 44, 2241 (1991); 

Review of the Homestake Solar Neutrino Experiment (^'^Cl Experiment): R. 
Davis, Prog. Part. Nucl Phys. 32, 13 (1994); 

For the Standard Solar Model, and various views on the status of the solar 
neutrino deficit, see, for example: J. N. Bahcall and M. M. Pinsonneault, Rev. 
Mod. Phys. 64, 885 (1992); 

S. Turck-Chieze et al., Phys. Rep. 230, 57 (1993); 

J. N. Bahcah, Nucl Phys. B (Proc. Suppl) 43, 41 (1995); 

S. Turck-Chieze and I. Lopes, Ap. J. 408, 347 (1993); 

T. J. Bowles and V. N. Gavrin, Annu. Rev. Nucl Part. Sci. 43, 117 (1994); 

S. T. Petcov, Nucl Phys. B (Proc. Suppl) 43, 12 (1995); 

K. V. L. Sarma, Int. J. Mod. Phys. A 10, 767 (1995); 

W. Kwong and S. P. Rosen, Mod. Phys. Lett. A 10, 1331 (1995); 

W. C. Haxton, Ann. Rev. Astron. Astrophys. 33, 459 (1995); 

and references therein. 

[24] G. Borner, The Early Universe: Facts and Fiction (Springer- Verlag, Berlin, 
1992). 

[25] R. Lieu et al., Science 274, 1335 (1996); 
S. Bowyer et al.. Science 274, 1338 (1996). 

[26] S. Bowyer (private communication. Winter 1996). 

[27] G. Schilling, Science 274, 1305 (1996); reporting on "a pair of papers soon to 
appear in Astronomy and Astrophysics" by H. Arp and H.-D. Radecke. 

[28] D. Buskulic et al., Phys. Lett. B 349, 585 (1995). 



24 



[29] R. M. Barnett et al. (Particle Data Group), Phys. Rev. D 54, 1 (1996). 

[30] K. Hirata et al., Phys. Rev. Lett. 58, 1490 (1987); 
R. M. Bioiita, Phys. Rev. Lett. 58, 1494 (1987); 

M. Aglictta et al., in The Standard Model, The Supernova 1987a, J. Tran 
Thanh Van, ed.. Proceedings of the Leptonic Session of the Twenty-Second 
Rencontre de Moriond (Editions Prontieres, Prance, 1987); 
E. N. Alexeyev et al., in The Standard Model, The Supernova 1987a, J. Tran 
Thanh Van, ed.. Proceedings of the Leptonic Session of the Twenty-Second 
Rencontre de Moriond (Editions Prontieres, Prance, 1987). 

[31] M. Roos, In The Standard Model, The Supernova 1987a, J. Tran Thanh Van 
ed., Proceedings of the Leptonic Session of the Twenty-Second Rencontre de 
Moriond (Editions Prontieres, Prance, 1993); 
J. N. Bahcah, and S. L. Glashow, Nature 326, 476 (1987). 

[32] E. Abies et al.. The MINOS collaboration, "P-875: A Long-baseline Neutrino 
Oscillation Experiment at Permilab," Pebruary 1995. 

[33] G. G. Padmabandu et al., Phys. Rev. Lett. 76, 2053 (1996). 

[34] E. P. Wigner, Ann. of Math. 40, 149 (1939). 

[35] V. V. Dvoeglazov, Int. J. Theor. Phys. 34 (1995) 2467. 
V. V. Dvoeglazov, Nuovo Cim. 108A (1995) 1467. 

[36] D. V. Ahluwalia, P, C, and T Structure of Space-Time (Kluwer Academic, 
Porthcoming). 

[37] D. V. Ahluwaha, M. B. Johnson, and T. Goldman, Phys. Lett. B 316, 102 
(1993). 



25 



